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Abstract 

We present OPE predictions for the hadronic mass and q 2 moments in inclusive 
semileptonic charmless B decays, with a lower cut on the charged lepton energy 
and an upper cut on the hadronic invariant mass. We include non-perturbative 
corrections through 0(l/mf) and perturbative contributions through O(a 2 (3o). We 
also investigate the range of the cut on the hadronic mass for which the local OPE 
can be considered valid and give estimates of the residual theoretical uncertainty. 

The study of these moments is important to constrain the effect of the Weak An- 
nihilation (WA) contributions and the shape of the distribution function, providing 
a more precise inclusive determination of \V u b\. 
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1 Introduction and General Motivations 



The determination of the V u b element of the CKM matrix, together with a sensible es- 
timate of the theoretical uncertainties involved in its calculation, is an important but 
challenging task. The magnitude of V u b, combined with V c b, determines the length of 
the left side of the Unitary Triangle, allowing for a stringent test of the Standard Model 
and a providing a powerful tool to detect hints of new physics pQ. While V c b, thanks to 
recent developments, has reached a theoretical uncertainty below 2%, the situation of V u b 
is much worse, with uncertainties about 10% 0]- Therefore improvements are needed 
here. 

The difficulties related with the determination of the V u b from inclusive semileptonic 
decays are well known. Inclusive decay rates can be calculated using a double expansion 
in a s (parton model) and AQCD/^b (Heavy Quark Expansion): 
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The measurement of B — > X u iv requires experimental cuts to suppress the large B — > 
X c iv background. Unfortunately, the introduction of kinematic cuts restricts the phase 
space available for the B — > X u iv decay, spoiling the global properties of the OPE. 

Several possibilities have been explored in the literature j3], involving cuts in the 
hadronic invariant mass (Mx), lepton energy (E{), leptonic invariant mass (g 2 ), light- 
cone component of the hadronic four-momentum (P+), and combinations of them. Each 
of these cases requires different techniques and a careful estimates of the uncertainties 
involved. 

In this talk, I will present a calculation of the moments in the hadronic invariant mass 
of the B — > X u iv distribution [5], and try to suggest the reason for its usefulness. Of 
course, as soon as the measurements of such moments will be available, one should verify 
their consistency with moments of B — > X c iv and B —>■ X^ in the OPE framework. 
The present combined fit jH] of all moments of inclusive B decays (hadronic mass, lepton 
energy) shows very good consistency of the available data and allows for the determination 
of the b quark mass and expectation values of the dominant power suppressed operators. 

However the main motivations for our analysis comes from recent experimental devel- 
opments. The high statistics accumulated at the B-factories by BaBar and Belle allows 
for the measurement, based on fully reconstructed events, of hadronic invariant mass 
distributions in inclusive B — > X u iv. This new generation of analysis can discriminate 
between charmless events and charmed background, even for relatively high values of 
hadronic invariant mass Mx- For instance, BaBar has been able to measure invariant 
mass distribution and first two moments with promising accuracy j7j . However, the exper- 
imental error is much smaller if one applies a cut on the hadronic mass M™* close to the 



kinematic boundary for charm production. On the other hand, theoretical calculations 
of inclusive B decays rely on the OPE, whose convergence is spoiled by severe cuts. A 
possible solution could be to raise M™* just enough to suppress non-perturbative effects 
that cannot be accounted for by the OPE (B meson distribution function effects). I will 
show a simple example in the Discussion (Sect. 4), but before let me briefly describe the 
content of our calculation (Sect. 2) and succinctly discuss the theoretical uncertainties 
involved (Sect. 3). 



2 Hadronic mass moments 

We determine the normalized integer moments of the squared invariant mass 

2 / dM\ M| n dY/dM 2 x 
{ x } ~ JdM 2 x dT/dM 2 x (2) 

and in particular the central moments defined as 

U X = (M X ) , U 2i3 = ({M 2 X - (M 2 X )) 2 ' 3 } . 

The physical hadronic invariant mass is related to parton level quantities by 

M\ = A 2 + 2m b AE + m 2 b s , 

therefore we can express the moments of M\ in term of combinations of moments of 
parton energy E Q and invariant mass of the parton s - Using the same approach, we can 
easily obtain also the moments of the distribution in q 2 . 

The building blocks that allow to construct both M\ and q 2 moments are: 

1 f d 3 T 

- t _i_ ^ r M^ G r _i_ p D n 4_ p ^f , a ^(i) i a2 s A) A (2) 

- %i) + + ^C(i,J) + ^D (iJ) + ^E (iJ) + -A {tj) + —A {hf) . 

This expression does not include non-perturbative terms of 0(l/mf), 0(a s /m 2 ), and 
perturbative corrections of 0(a 2 ). 

The moments of Eq and sq and of their product are obtained in the local OPE and 
are expressed in terms of the heavy quark parameters. 

We employ the "kinetic" scheme jHj (Wilsonian scheme with a hard factorization scale 
fi ~ 1 GeV). We start from on-shell expressions, and express the on-shell parameters in 
terms of the /z-dependent "kinetic" parameters: kinetic mass of the b quark m^/i), kinetic 



expectation value p 2 (p), and Darwin expectation value p%(p)- At p — >0 one recovers the 
results of the on-shell scheme. 

In the case of charmless decay, we introduce an additional parameter X M = 8 In ml / p\ . 
It is well known [TU| that the coefficient function -D(o,o) ( see Eq. EJ) has a logarithmic 
divergence as m q — > 0, i.e. contains a term proportional to In (ml/m 2 ). This problem 
is solved when we include the contribution of the Weak Annihilation (WA) operator. In 
fact, a one- loop penguin diagram that mixes the four-quark operator into p z D replaces 
lum^/ml by In p\ q /mi, where p^ q is the normalization point of the WA operator. Vary- 
ing the parameter X M allows us to estimate the effects of WA. The associated variation in 
the results for the moments accounts for the non-valence piece in the expectation value 
of the WA operator (flavor singlet WA). This contribution does not distinguish between 
B + and B° (non-singlet WA effect). 

The calculation is implemented in a FORTRAN code 1 and includes non-p ert urbat ive 
corrections 0(1/ ml) [TT] and 0(1/ ml) 9J, and perturbative corrections 0(a s ) and 
0(a 2 s Pq) [12]. Partial checks have been performed to reproduce previous results [T3j . 

3 Theoretical Uncertainties 

There are several sources of theoretical uncertainty in our predictions for the moments. 
Let me briefly list them and describe the methods used to estimate their size. We should 
consider: (i) Uncalculated 0(a 2 s ) and perturbative corrections to the Wilson coefficients; 
(ii) Missing 0(1/ m^) non-perturbative effects; (iii) Error from the scale in X M ; (iv) WA 
contributions; (v) Effects of Fermi Motion. 

We estimate the effect of missing higher order corrections (i) and (ii), by varying 
the parameters p 2 , p 2 Gl a s , p 3 D , p\g, and m& in an uncorrelated way. Alternatively, we 
can estimate the size of missing higher orders comparing the results obtained for the 
moments using different rearrangements of non-perturbative corrections. As I mentioned 
previously, the effects of (iii) and (iv) are closely related [3]: varying p± q over a reason- 
able interval, we estimate the effect of the flavor singlet WA. Those contributions are 
concentrated at maximal q 2 (small hadronic invariant mass) and are therefore suppressed 
in the moments of M\. Finally, we estimate the effect of Fermi motion (v), introduced 
by the Mx cut, by smearing the tree- level differential rate with an exponential distri- 
bution function. This estimate depends on the functional form adopted, for instance an 
exponential form leads to more conservative estimates than a Gaussian ansatz. While 
find relatively small effects for the first moment, the second and the third one are more 
sensitive. A more complete and satisfactory treatment of these effects is required, and it 
is currently under investigation. 

1 Available upon request from the authors of Ref. [5]. 



4 Discussion and Conclusions 



As an example of the possibilities offered by the analysis based on the moments, we 
present an estimate of the error on a the determination of mj based on U\ only. We con- 
sider a mild M\ <5.6 GeV cut, that has the advantage of removing the experimentally 
most poorly known region of phase-space for charmless decays, reducing the experimental 
error by almost a factor two, respect to the case without the cut The experimental 
uncertainty, considering the statistics that can be accumulated in the near future, can 
be estimated as 5U^ P = 11% [14 . On the theoretical side, considering the sources of 
uncertainty listed above, we obtain dUf 1 = 9% at M™* =2.5 GeV[5 . The combined ef- 
fect of theoretical and experimental uncertainties, together with the use of the linearized 
formulas of Ref . |x , roughly leads to ~ 80 MeV. 

To summarize, we presented predictions for the first three moments of the hadronic 
invariant mass distribution in B — > X u £u, that includes corrections through O(a 2 s fio) 
and 0(1 /m\) and cuts on the lepton energy and on the invariant hadronic mass. The 
calculation is implemented in a FORTRAN code, available upon request. The work is 
in progress: we are currently improving our code, allowing all possible combined cuts 
(including q 2 and light-cone variables). A lot is still to be done, in order to constrain 
WA (singlet/non-singlet) and better understand Fermi motion effects. An analysis based 
on moments with high M™* is challenging for experiments, but reduces the theoretical 
uncertainties related with OPE predictions. The preliminary results are promising. 
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